Quantum blockade and loop current induced by a single lattice defect in graphene 

nanoribbons 
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We investigate theoretically the electronic transport properties in narrow graphene ribbons with 
an adatom- induced defect. It is found that the lowest conductance step of a metallic graphene 
nanoribbon may develop a dip even down to zero at certain values of the Fermi energy due to the 
defect. Accompanying the occurrence of the conductance dip, a loop current develops around the 
defect. We show how the properties of the conductance dip depend on the parameters of the defect, 
such as the relative position and severity of the defect as well as the width and edges of the graphene 
ribbons. In particular, for metallic armchair-edges graphene nanoribbons, whether the conductance 
dip appears or not, they can be controlled by choosing the position of the single defect. 
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I. INTRODUCTION 



Graphene has attracted intensive interest recently 
because of its novel properties, such as the room- 
temperature quantum Hall effect [l| and the Dirac- 
equations-governed electrons Q 01 . Moreover the proper- 
ties of the graphene nanoribbons [5[ change notably due to 
various edges p, different symmetries [J, constrictions 
and even controllable defects by irradiating the material 
with electrons and ions @ [HI • Any wanted structures 
could be achieved by engineering the graphene with the 
help of state-of-art fabrication technologies, which sup- 
ports the graphene as a prospective candidate for signif- 
icant applications in electronics- and magnetics-related 
devices. 

Defects have been revealed to be an influential fac- 
tor in determining the transport properties of structures 
in nanoscale. For example, it has been known that de- 
fects can reduce or totally block the quantum conduc- 
tance in carbon nanotube [i2|[3[3- Graphene could 
be regarded as an unrolled carbon nanotube, and hence 
the transport of graphene in the presence of defects is 
also an interesting problem. Actually, several kinds of 
defects in a single layer of graphene have already been 
paid attention to, such as pe ntagon / heptagon defects 
im [iGj , vacancies [13] , adatoms jl8l| |T9j , substitution 2C | , 
disorder [2lj and even combinations of some of them 22 1. 

Although the electronic energy band structures and 
magnetic properties in the presence of defect have been 
investigated in both graphene and graphene ribbons, how 
these properties manifest in transport is still lacking in 
attention. Recently, it has been found that the interior 
pentagon/heptagon defects or the defects such as va- 
cancies or disorder at the edges [ll] would lead to total 
quantum blockade in graphene nanoribbons at certain 
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Fermi energies. However, how the quantum blockade is 
connected with the severity or the positions (not limited 
to the edges) of the defect, and the size or the symme- 
try of the graphene nanoribbons is not so clear. As the 
adatom defect would not change the topological struc- 
ture of graphene lattice to some extent, the connection 
between these factors and the defect can be unambigu- 
ously understood by comparing with the clean graphene 
system. As we will show, all these factors are really crit- 
ical on the influence of the defect. Therefore, a careful 
study of the influence of the adatom defect on the trans- 
port properties of graphene nanoribbons is deserved. 

In this paper, we investigate the electronic transport 
in narrow graphene ribbons in the presence of a single 
defect, which could be induced by the absorption of a hy- 
drogen atom. The model of the defect is based on the ab- 
initio calculation. We find that the conductance profile 
is greatly reshaped by the defect. Compared with the in- 
teger steps of perfect graphene nanoribbons, the conduc- 
tance in the presence of a single defect drops even down 
to zero in the low energy regime, leading to a transfor- 
mation from metal to semiconductor. Accompanied by 
this global conductance blockade, a loop current builds 
up around the defect. The conductances of the defec- 
tive graphene nanoribbons with different widths, edges, 
and positions of the defect are studied. We find the con- 
ductance blockade and the coexistence of the loop cur- 
rent sensitively depend on the severity and relative posi- 
tion of the defect. For metallic armchair-edges graphene 
nanoribbons, whether the conductance blockade and the 
coexistence of the loop current develop or not, they could 
be controlled by the position of the defect. 



II. MODEL AND RESULTS 

The transport properties in graphene ribbons under 
low bias voltages are governed by the tt bonded elec- 
trons, which are modeled by the tight-binding spinless 
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Hamiltonian 



H 



Cj + h.c. 



(1) 



where c| (ci) creates (annihilates) an electron on the site 
i, t is the nearest-neighbor hopping energy and is set as 
2.7 eV in this paper. The distance between the nearest 
neighbor carbon atoms is a (~1.42 A). 

The conductance can be calculated by the Landauer- 
Buttiker formula based on the Green function 
methodfiilfiijllilfi^]. The formalism separates 
the entire system into three parts, the central re- 
To mimic the transport through 
two semi- 



entire 
gion and two leads. 

infinitely long ribbon with specific edge, 
infinite graphene leads are connected to the central 
region, see FigH] (a). The zero temperature con- 
ductance G of the graphene nanoribbons is given 
by G{E) = ^TT[TsiE)G''{E)TD{E)G''{E)], where 

G^(-\E) = {E I-n- - S^^'V is the re- 

tarded (advanced) Green-function matrix of the central 

part, (i?) is the retarded (advanced) self-energy 

due to the source (drain) [labeled as S(D) in the Fig[T] 



(a)], and Tsm)iE) 



^S(r>)(^) 



Here I 



represents the unitary matrix. The local current density 
at the Fermi level E between two neighboring sites i and 
j can be expressed as 



tjG^iiE)], 



(2) 



where G" — GTsG" is the electron correlation func- 
tion and Hij is the corresponding matrix element of the 
Hamiltonian. 

In the absence of the defect, the well-known conduc- 
tance of the perfect graphene nanoribbon with zigzag 
edges is shown by the dot line in Fig. 1(b). Due to 
the metallic nature of the system, there exists no energy 
gap. The quantum conductance begins with Go (defined 
as 2e^/h) and displays a typical step-increasing feature 
as a function of the Fermi energy. We also show in Fig. [T] 
(a) the spatial distribution of the local current density at 
a selective energy value of E=0.32t. One can see that the 
current throughout the nanoribbon is uniform along the 
transport direction and forms unambiguous paths. At 
this energy, the current density amplitude at the central 
region of the nanoribbon is bigger than that at the edge 
regions. 

Now, let us focus on the effect of a single defect site 
on the quantum conductance and local electric current 
distribution in the graphene nanoribbons. To physically 
illustrate how to model the presence of a single defect site 
by coupling a single carbon atom in the graphene with an 
adatom atop, we have simulated the atomic hydrogen ad- 
sorption on the graphene [see Fig. [DJa)] at a low coverage 
of 0.02 monolayer (namely, one hydrogen atom presents 
in a 5x5 graphene supercell), by using density- functional 
theory within the local density approximation and a su- 
percell approach. The total energy calculation shows that 




FIG. 1; (Color online) (a) The local electric current distribu- 
tion in the perfect graphene ribbons. The sizes of the arrow 
represent the magnitude of the electric current between any 
two neighboring sites, which are linearly normalized to the 
maximum value. The source and drain are labeled with S and 
D respectively, (b) The conductances of a perfect graphene 
nanoribbon (dot line) and the one with a defect (solid line). 
The position of the defect is shown in (c). (c) The local 
electric current distribution in the graphene nanoribbon with 
a defect marked by a hollow circle. The Fermi energy is 
E = 0.32t. 



the most stable adsorption site for H adatom is the on-top 
site and the H-C bond is characterized by a hybridization 
of H Is and C 2pz states. Prominently, it is found that 
the C atom beneath the H adatom is pulled a little bit out 
of the garphene plane by H-C chemical bonding, which 
obviously alters the coupling between this C atom and its 
three neighboring C atoms and thus results in a single lat- 
tice defect site. For further clarification, Fig. [DJb) plots 
the calculated electron density difference An(r), which is 
obtained by subtracting the electron densities of nonin- 
teracting component systems, rigraphonc (r) -fnH(r), from 
the density n{r) of the H/graphene system, while retain- 
ing the atomic positions of the component systems at the 
same location as in the H/graphene system. One can see 
that the charge redistribution only occurs around the H- 
C bond and influences the three neighboring C-C bonds. 
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(b) 



FIG. 2: (Color online) Adsorption of H on graphene: (a) The 
side view of H/ graphene system with an H adatom (the small 
green ball) adsorbed at on-top site of the graphene (the big 
gray balls), (b) The top view of the isosurface of the electron 
density difference (±0.05e/A'^). Charge flows from the blue 
into the red regions. 
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FIG. 3: (Color online) Left: Conductances with the ratio td/t 
decreasing from top to bottom. The system is the same one as 
in Fig[T] Right: Conductances of the graphene nanoribbons 
with different widths. In this case, the position of the defect 
and length of the ribbon is the same as the one in Fig[T] 
The width changes by cutting or adding chains of C atoms 
at the bottom of the graphene nanoribbons. L is number of 
the atoms connecting the source lead. For every plot, the 
conductance of the perfect graphene nanoribbon is shown for 
reference with the dotted lines. 



while the charge densities of other C atoms are not af- 
fected at all. 

Based on the physical picture revealed in Fig. [21 we 
model the presence of a single defect lattice site by mod- 
ifying the interactions between this lattice site and its 
nearest neighboring sites while keeping other hopping el- 
ements unchanged. Therefore we change the hopping 
amplitude between the single defect site [hollow circle 
in Fig. [Ijc)] and its three nearest neighbors from t to, 
for example, td—Q-it. The conductance of the system 
with the defect is calculated and plotted with the solid 
line in Fig. [Hb). One can see that the conductance is 
prominently reshaped by the presence of the single de- 
fect site. The explicit changes include: (i) compared to 
the perfect case, the quantum step feature disappears; 
(ii) the conductance drops for most Fermi energies, (iii) 
zero-conductance dips appear at certain Fermi energies. 
These defect-induced phenomena are understandable be- 
cause the defect is a kind of destructive factor for the 
electron transmission in otherwise perfect graphene lat- 
tices. We now concentrate on the conductance dip at 
E=0.?>2t. Since the conductance drops to zero at this 
point, it means that the transmission is totally block- 
aded by the defect. The dip also implies the appear- 
ance of energy gap and consequently a transformation 
from metal to semiconductor happens. For further ex- 
planation to the conductance dip, we calculate the local 
current density at the Fermi energy (i?=0.32t) and plot 
the result in Fig. [T](c). Compared to the ideal graphene 
nanoribbon system [Fig. [Ija), plotted at the same energy 



i?=0.32t], it shows in Fig[T] (c) that the pattern of local 
current density distribution is prominently changed in 
the defect-included system and there appears loop cur- 
rent circling around the defect site. The magnitude of 
this loop current around the defect can be larger than 
the maximal current in FiglTJa). Note that the current 
is normalized to the maximum value for each case. The 
loop current would produce an effective magnetic field, 
which is consistent with t he p revious experimental and 
theoretical predication. [rTj [l8| [28t 

The formation of the loop current could be explained as 
the interference of all transport channels. In the two di- 
mensional graphene lattices, an electron could jump from 
source to drain via many different transport channels. 
The local electric current actually demonstrates the in- 
terference results of these channels. A defect would cause 
the blockade of relative transport channels to some extent 
and therefore change the local electric current. A discus- 
sion using a simple one-dimensional model is given in 
reference .23j. As the different transport channels are nu- 
merous in the two dimensional graphene lattices, whether 
the loop current appears and what kind of patterns it 
takes depend on the whole environment of the defect, 
such as the edges and width of the graphene ribbons, the 
position and severity of the defect, etc. Therefore, we 
would dwell on the influence of these parameters on the 
transport properties in the defective graphene ribbons. 

We can define the severity of the defect site by the ra- 
tio td/t. The vacancy defect is actually the defect with 
td/t=0. A question thereby arises: How effective could 
the defect cause the conductance to drop to zero? In 
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FIG. 4; (Color online) Conductance calculated for the case 
when the on site energy of defect adds 0.1 t, shown with 
solid line. Other parameters are the same as Fig[ljb). For 
reference, the case of a perfect lattice (dot line) and the one 
with a single defect but assuming on site energy unchanged 
(dashed line) are also shown. 



the following, we will show that it depends on the com- 
petition between the width of the graphene ribbon and 
the severity of the defect. For a narrow graphene rib- 
bon, the zero conductance exists in most cases when for 
small ratios of td/t. In the left panel of Fig[3l we present 
the conductance of the same defective graphene ribbon 
as in FigH] but with different td/t. We can see that as 
the ratio td/t decreases, the dip begins to appear in the 
lowest conductance step and conductance reaches zero 
at td/t=QA. This critical value at which zero conduc- 
tance is obtained is also dependent on the width of the 
graphene ribbon. To show the influence of the width, we 
modify the width of the defective graphene nanoribbon 
shown in FiglT] by cutting or adding chains of C atoms 
on the lower edge of the ribbon. The width L is denoted 
by the number of chains. In the right panel of FiglSl the 
conductances of the graphene ribbons with four different 
widths are plotted. It is obvious that as the width of 
the graphene ribbon gets wider, the dips become shallow 
and narrow and finally almost disappear. For the case 
of L = 38, the values of conductances almost get equal 
between E = Q.22t and, for example, E — Q.lt. However, 
the loop current only exists near E — 0.22i, not E ~ O.lt, 
although the net current cannot be totally blocked by the 
defect. 

Besides the change of inter-site hopping matrix ele- 
ments by the coupling of a single lattice site with the 
adsorbed atom, the on-site energy of this lattice site may 
also be changed due to the bonding and anti-bonding 
hybridizations. The conductance by changing the on-site 
energy of the defect site by Q.lt is shown in FigUl We 
see that the conductance dip still appears but lies at a 
shifted energy not far away. The distance depends on 
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FIG. 5: (Color online) (a) The graphene nanoribbon with 
four different defect sites to be considered respectively, which 
are labeled with A, B, C and D. (b) Conductances of the 
graphene nanoribbon with one defect located at one of these 
four positions respectively. Inset is the amplification of the 
part where the conductances drop to zero. 



the changed amplitude of on-site energy. At the conduc- 
tance dip, the loop current also exists without apparent 
difference from the one shown in Fig[l] (c). So we con- 
clude that the on-site energy is not an influential factor 
to the occurrence of conductance dip and formation of 
loop current induced by the single defect in this model. 

We have also studied the influence of the position of the 
defect site on the conductance of the graphene nanorib- 
bons. It is found that the horizontal shift of a single de- 
fect produces no changes to the conductance, while the 
different defect positions in vertical direction leads to dif- 
ferent conductances and hence the center of loop current 
patterns. In Fig. we choose respectively four sites 

(labeled by A, B, C and D) as a single defect site to make 
a comparison. The conductance corresponding to these 
four choices for the defect site are plotted in Fig. [^b). 
One can see that for each case the conductance drops to 
zero with a dip. The width of the dip becomes smaller 
with the site going closer to the boundary, see the inset 
in FigO (b). By comparing the current flowing through 
these sites in perfect lattice, we see that the defect lo- 
cating on the position with bigger current in the perfect 
lattices would cause a conductance dip with a broader 
width. 

The role of the presence of a single defect lattice site be- 
comes more important in the graphene nanoribbons with 
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FIG. 6: (Color online) (a) Local electric current distribution 
(arrows) in a metallic armchair-edges graphene ribbon. Sites 
A (left-triangle) and B (right-triangle) represent two kind of 
defect to be considered respectively, (b) Conductance of the 
system with the defect located at position A. (c) Conductance 
of the system with the defect located at position B. The con- 
ductance of the system without any defect is also given with 
dotted lines for reference in (b) and (c). 

armchair edges. As we know, a graphene nanoribbon 
with armchair edges has the zero energy gap only when 
the number of atoms on a zigzag edge connecting the 
source (or lead) is 3M— 1 (M is an integer). This kind of 
graphene nanoribbons hence has a quantum conductance 
beginning from Go with increasing the electron Fermi en- 
ergy. We plot the local current distribution of the system 
without the defect in Fig. [HJa), which clearly displays a 
series of well-separated pathes for current flow. Each 
path has the same weight in current amplitude. It is the 
very result of interference between the transport chan- 
nels in the metallic armchair-edges graphene ribbons. [1^ 
According to this stripe current distribution, the lattice 
sites in the graphene nanoribbon can be therefore classi- 
fied into two parts, i.e., sites like A [see Fig. [6]^ a)] carry- 
ing current flow and sites like B without carrying current 
flow. We set the single defect site to be on site A or 
on site B, and then calculate the corresponding quantum 
conductance respectively. The results are shown in Fig. 



[Ijb) and ^c) , respectively. Remarkably, it is shown in 
Fig. [6jb) that when a single defect is placed on site A, 
no conductance dip occurs at all and the first quantum 
conductance step is invariant upon the presence of the 
defect. Whereas, when the defect is placed on site B, 
two dips develop in the conductance spectrum, as shown 
in Fig. [Gjc). Again, at the energy position E—Q.32t of 
the conductance dip, the spatial distribution of the inter- 
nal current (not shown here) is featured by a loop current 
around the defect site. Thus, the lattice sites B are more 
important than the lattice sites A for the transport in 
this kind of armchair-edges graphene ribbons. This is an 
important ingredient for practical control of the graphene 
transport by manipulating the position of a single defect 
site. From the point of the interference among different 
transport channels mentioned above, the defect in site B 
gives a kind of destructive interference. The local cur- 
rent density distribution pattern of the perfect metallic 
armchair-graphene ribbons indicates that the total con- 
ductance is mainly contributed by several parallel and 
independent channels. At the site A, the local density is 
zero. In other words, the site A is decoupled and isolated 
from these channels. That is why the graphene ribbon is 
immune to the defect on site A. 



III. CONCLUSION 

In summary, we have studied the influence of a sin- 
gle defect on the electron transport properties in narrow 
graphene nanoribbons with different edges, widths and 
positions of the defect. It has been found that the con- 
ductances may drop even down to zero at certain Fermi 
energies in the presence of a single defect site for metallic 
graphene nanoribbons. The conductance dip is accom- 
panied by the formation of a remarkable loop current 
around the defect. The width of the conductance dip 
is sensitive to the position as well as to the severity of 
the defect site. For the graphene nanoribbons with arm- 
chair edges, it has been found that the transport prop- 
erties, including the low-energy conductance step, the 
conductance dip and the loop current, can be controlled 
by choosing the position of the single defect. This re- 
sult enlightens us on controlling the conductance of the 
graphene nanoribbons by simply manipulating the de- 
fects due to the absorbed atom. 
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